Dynamical creation of bosonic Cooper-like pairs 
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We propose a scheme to create a metastable state of paired bosonic atoms in an optical lattice. 
The most salient features of this state are that the wavefunction of each pair is a Bell state and that 
the pair size spans half the lattice, similar to fermionic Cooper pairs. This mesoscopic state can be 
created with a dynamical process that involves crossing a quantum phase transition and which is 
supported by the symmetries of the physical system. We characterize the final state by means of a 
measurable two-particle correlator that detects both the presence of the pairs and their size. 



Pairing is a central concept in many-body physics. It 
consists on the existence of quantum or classical correla- 
tions between pairs of components of a many-body sys- 
tem. The most relevant example of pairing is BCS su- 
perconductivity. In a superconductor, an attractive in- 
teraction causes electrons to organize into Cooper pairs, 
bosonic quantum objects where electrons are perfectly 
anticorrclatcd in momentum and in spin. In the language 
of second quantization this is described by the BCS vari- 
ational wavefunction 



\<hcs) =H(u k + v k A{)\0), 
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where A k = c^cL^ * s an °P cra t° r that creates one such 
Cooper pair. Remarkably, the fact that pairing occurs in 
momentum space means that the constituents of the pairs 
are delocalized and share some long-range correlation. 

Nowadays, pairing is also a key research topic in the 
field of ultracold atoms. While atomic interactions are 
short range, they can be enhanced using Feschbach reso- 
nances. These resonances allow both to create Cooper 
pairs of fermionic atoms [3, 0, 0| and to observe the 
crossover from these large, delocalized objects to a con- 
densate of bound molecular states. Realizing similar ex- 
periments with bosons is difficult, because attractive in- 
teractions can induce collapse. In the paper by Paredes 
[H it is suggested to load an optical lattice with hard- 
core bosonic atoms in two internal states that attract 
each other so as to induce BCS-like pairing. A different 
approach, followed in Ref. Q, is to load an optical lattice 
with pairs of atoms in the regime of strong repulsive in- 
teractions. Such metastable on-site pairs are robust and 
survive about 700 ms in the lattice. 

In this Letter we propose a method to create long-range 
paired states of bosonic atoms using entangled states as a 
resource. The method uses an optical lattice of arbitrary 
geometry which is loaded with entangled bosons in an 
insulator state. We will consider both on-site pairs 
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such as the ones demonstrated in Ref. [6[ and a larger 
family of singlet and triplet states 



L/2 

II 4 , 



JO), A, 



(3) 



some of which have been created using optical super- 
lattices 0, H| ■ We propose to dynamically increase the 
mobility of the atoms, entering the superfluid regime. 
During this process, the pairs will grow in size until they 
form a stable gas of long-range Cooper-like pairs that 
span about half the lattice size. Contrary to works on the 
creation of squeezed states 0] , the evolution considered 
here is not adiabatic and the survival of entanglement is 
ensured by a symmetry of the interactions. 

This paper is organized as follows. First, we present 
the Hamiltonian for bosonic atoms which are trapped in 
a deep optical lattice, have two degenerate internal states 
and spin independent interactions. Next, we prove that 
by lowering the optical lattice and moving into the super- 
fluid regime, the Mott-Bcll entangled states (0)-([3]) evolve 
into a superfluid of pairs. We then introduce two corre- 
lators that detect the singlet and triplet pairs and their 
approximate size. These correlators are used to interpret 
quasi-exact numerical simulations of the evolution of two 
paired states as they enter the superfluid regime. Finally, 
we suggest two procedures to measure these correlators 
and elaborate on other experimental considerations. 

We will study an optical lattice that contains bosonic 
atoms in two different hypcrfine states (a =t , J.)- In the 
limit of strong confinement, the dynamics of the atoms 
is described by a Bose-Hubbard model 
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Atoms move on a d-dimensional lattice (d = 1, 2, 3) jump- 
ing between neighboring sites with tunneling amplitude 
J CT , and interacting on-site with strength U aa i . The Bosc- 
Hubbard model has two limiting regimes. If the interac- 
tions are weak, U <C J, atoms can move freely through 
the lattice and form a superfluid. If interactions are 
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strong and repulsive, U ^> J, the ground state is a Mott 
insulator with particles pinned to different lattice sites. 

As mentioned in the introduction, we want to design 
a protocol that begins with a localized entangled state 
©-([2]), and, by increasing the mobility of atoms, enlarge 
these pairs until the final state describes a gas of general- 
ized Cooper pairs of bosons. In our proposal we restrict 
ourselves to symmetric interaction and hopping ampli- 
tudes 

U = U n = U U = U n > 0; J = J T = J i > 0. (5) 

This symmetry makes the system robust so that, even 
though bosons do not stay in their ground state, they 
remain a coherent aggregate of pairs, unaffected by colli- 
sional dcphasing. We will formulate this more precisely. 

Let us take an initial state of the form given by either 
Eq. (0) or If we evolve this state under the Hamilto- 
nian with time- dependent but symmetric interaction 
U aa i = U(t) and hopping J a — J{t), the resulting state 
will have a paired structure at all times 

|^))-£ C (t;k)4 ifc2 ...4 2iifc2 j0), (6) 
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where c(i, k) are complex coefficients to be determined. 

The proof of this result begins with the introduction of 
a set of operators Cy := 2 CT c L c j<r which form a simple 
Lie algebra [Cy,Cfc;] = GuSjk — CkjSu. The evolution 
preserves the commutation relations and maps the group 
onto itself. This is evident if we rewrite the Hamiltonian 

(id) * 

The evolution operator satisfies a Schrodinger equation 
ih^V(t) = H(t)V(t), with initial condition V(0) = I. 
Since the Hamiltonian only contains CV, operators we 
conclude that V(t) is an analytic function of these gen- 
erators. Let us now focus on the evolution of state ([3]), 

given by \f(t)) = V(t)l\U 2 i A-i^M ■ We wil1 use the 
commutation relations between the generators of the evo- 
lution and the pair operators [Aij, Cm] = SikAij + SjkAu, 
which are valid for any of the pairs in Eq. ([3]). Formally, it 
is possible to expand the unitary operator V(t) in terms 
of the correlators Cjj and commute all these operators 
to the right of the A's, where we use CV,|0) = 0. After 
doing this one is left with Eq. ([6]). A similar proof applies 
for the on-site pairs ((2]) and for any initial state which 
already has a paired structure ([6|). 

The previous result includes a very simple case, which 
is the abrupt jump into the non-interacting regime, U = 
0. This problem is integrable and for the initial condi- 
tions © and ^ the evolved state can be written as 

N 
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The wavepackcts w(i,j,t) form an orthogonal set of 
states, initially localized w(i,j,0) oc <5y or w(i, j, 0) oc 
Sij+i and for large times close to a Bessel function [181 ] . 
We remark that the pair wavcfunctions ((6]) and (|8|) can 
include valence bond states, however represent a gener- 
alization both in the fact that pairs can overlap or be 
triplets. 

In a general case, computing the many-body pair wave- 
function, c(t; k), is an open problem. Nevertheless we can 
prove that the final state does not remain in the ground 
state in the superfluid regime, no matter how slowly one 
changes the hopping and interaction. For the states in 
((3|) this is evident from the lack of translational invari- 
ance. Let us thus focus on the state ([2]) generated by 
An = Ci\ , which has an equal number of spin-up and 
down particles N^i = N/2. The ground state of the 
same sector in the superfluid regime, U = 0, is a number 
squeezed, two-component condensate [§] 

iv>NN)«4f /2 4r /2 io>> (9) 

with coo- = -7= ^2f—i Ci a . Note that we can also write this 
ground state as an integral over condensates with atoms 
polarized along different directions 

|Vnn) oc J dee-^ic^ + e^c^flO), (10) 

When this state is evolved backwards in time, into the 
J = regime, each condensate transforms into an insu- 
lator with different polarization yielding 

ivw^> E n( c l T ) nfc (4i) 2 " nfc io). (11) 

n,J2n k =N/2 k 

Since this state is not generated by the An = CijCjj op- 
erators, we conclude that this particular state ([2]), when 
evolved into the superfluid, leaves the ground state. Fur- 
thermore, since different pairs in Eq. ([2]) are related by 
global rotations, this statement applies to all of them. 

For the rest of this Letter we focus on two important 
states: the triplet pairs generated on the same sitelgl and 
the singlet pairs generated on neighboring sites [7|, [8J , 

H't) = n|(cl T 2 + clf)|0>, and (12) 
i=l 
L/2 

\^s) = -^=(c 2l -iTC2ii — C 2l -14.C2it)|0>i ( 13 ) 
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respectively. Our goal is to study the evolution of these 
states as the mobility of the atoms is increased, suggest- 
ing experimental methods to detect and characterize the 
pair structure. The main tools in our analysis are the 
following two-particle connected correlators 

G% := (ct T ct T CijCjM^ - (c^cnj (c]t c ji) > ( 14 ) 
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combined in two different averages 
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and what we call the pair size 
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A variant of the correlator G T has been used as a pair- 
ing witness of fermions . We expect these correlators 
to give information about the pair size and distribution 
also in the superfluid regime. This can be justified rig- 
orously in the case of an abrupt jump into the super- 
fluid, where the pair wavepackets remain orthogonal and 
where Ga and R characterize the spread of the wavefunc- 
tions w(i,j,t). First, note that the single-particle expec- 
tation values such as (cj^Cij) are exactly zero since Nf 
and iVj, arc even for the triplet state ipr and balanced for 
the singlet state ips- Second, the two-particle correlators 
only have nonzero contributions where the destruction 
and creation operators cancelled and subsequently cre- 
ated the same pair. Combining Eqs. (TH|) and © gives 

G l = ^2\w(i-n,j-n,t)\H, (17) 

n 

G fj = -EW-v-M)l 2 , 

n 

where we have used the symmetry of the wavefunction, 
w(i,j,t) = w(j,i,t). Particularized to the initial states, 
the triplet ipr gives Gfj = Sij, G^ = <5ao> G t = 1, and 
R T = 0, as expected from on-site pairs. The singlet 
pairs described by tps, on the other hand, yield Gfj = 

I^Aij G s = — i, and initially 

1. 

For a realistic study of the evolved paired states we 
have simulated the evolution of tpT and ips under the 
Bose-Hubbard model as the hopping increases nonadia- 
batically in time 



— 2(^+1 g a — —5. 

span two sites R s 



J{t) = vx (tU/h) x U, 



(18) 



with ramp speeds v = 0.5, 1 and 2 in adimcnsional units. 
The simulations were performed using Matrix Product 
States (MPS) on one-dimensional lattices with up to 20 
sites and open boundary conditions 0, E3| ■ Given the 
small size of these lattices, we expect these simulations to 
appropriately describe even the superfluid regime, where 
the size of the energy gaps and the high occupation num- 
ber per site make the MPS simulation more difficult. 

For each of these simulations we have computed all cor- 
relators and pair sizes, values which are plotted in Fig.Q] 
Let us begin with the triplet pairs: initially the only rele- 
vant contribution is the short-range pair correlation Gq ; 
then the pair size increases monotonously from R = 




FIG. 1: We plot the (a) triplet and (b) singlet correlators 
for the evolution of ipr and ips respectively, in a lattice of 
L — 20 sites. The circles, triangles and squares denote G, Go 
and their difference, (c) Pair size R for the singlet (line) and 
triplet (cross) states, for a ramp speed v = 0.5, 1 and 2 (solid, 
dash, dash-dot). The vertical line J/U — 1/3.84 marks the 
location of the phase transition. (Color online) 



up to R ~ L/2, where it saturates. At this point, the 
pairs have become as large as the lattice permits, given 
that we have uniform density and open boundary con- 
ditions. The singlets have a slightly different dynamics. 
The antisymmetry of the spin wavefunction prevents two 
bosons of one pair to share the same site and thus R = 1 
initially. However, this antisymmetry seems also to affect 
the overlap between pairs, as it is evidenced both in the 
slower growth R(t) and in the smallncss of Gq . 

Concerning the speed of the process, we have simu- 
lated ramps over a timescale which is comparable or even 
shorter than the typical interaction time, 1/U, so that the 
process is definitely not adiabatic. Nevertheless, the pairs 
seem to have enough time to spread over these small lat- 
tices. Note also that the spreading of atoms begins right 
after the value J/U ~ 1/3.84 where the one-dimensional 
Insulator-Supcrfluid phase transition takes place 

The system of delocalized Cooper-like pairs can also be 
regarded as a mean of distributing entanglement in the 
optical lattice. Following this line of thought we have 
used the von Neumann entropy to measure the entangle- 
ment between two halves of the optical lattice. Apart 
from a logarithmic contribution induced by the dclocal- 
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ization of particles, which is also present in a single- 
component condensate, there is an additional contribu- 
tion caused by the spreading of pairs across the lattice. 
However, this contribution does not reach a magnitude 
0(N/2) which corresponds to perfectly splitting N pairs 
between both lattice halves. We conjecture this is due to 
the pairs being composed of distinguishable particles. 

The pairing correlators G T ' can be decomposed into 
density-density correlations and measured via the noise 
interferometry technique developed in Ref. [IH and ap- 
plied in Ref. To prove this, let us introduce the 
Schwingcr representation of angular momenta 



S x (i) 

Sy(i) 

S.(i) 



71 ( c h c ii 
i(4 T c lT - cj^cij.) = i(n lT - Tin). 



(19) 



We will focus on the real part of the correlation matrices, 
defined as Gij = 2Re(G,j). For the family of states con- 
sidered here this matrix is the same as the original one 
[See Eq. (fl~7|) ] and it is related to simple spin correlations 

±(S x (i)S x (j))T^(S y (i)S y (j)), (20) 

±(S,(Z)) (S X (j)}±±(Sy(i)) (Sy(j)). 
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We now introduce two global rotations in the hyper- 
fine space of the atoms U X:V = exp [±«-f ^2 k S ViX (k)~\ . 
The operators U x and U y take the S x and S y opera- 
tors into the S z , respectively. These rotations can be 
implemented experimentally without individual address- 
ing. Using these unitaries we can relate the previous spin 
correlations to simple density operators. For instance 

(S x (i)S x (j)) = \ (UUmi - nn){n j} - n n )U x ) , (21) 

shows that the S X S X arises from all possible density cor- 
relations after applying a n/2 pulse on the atoms. 

Another possibility is to apply the ideas put forward in 
Ref. [13] • These methods rely on the interaction between 
coherent light and the trapped atoms to map quantum 
fluctuations of the atomic spin onto the light that crosses 
the lattice. Using this technique it should be possible 
to measure both the single-particle and the two-particle 
expectation values that constitute G T S . 

Experimental imperfections are expected not to affect 
the nature of the final state. The influence of stray mag- 
netic and electric fields can be cancelled by working with 
the singlet pairs, which are insensitive to global rotations 
of the internal states and have large coherence times. 
More important could be the influence of any asymme- 
try in the interaction constants. However, assuming this 
asymmetry to be of the order of 1%, the effect can only 
be noticeable after a time t = lOOh/U, which is longer 
than the evolution times suggested here. 

Summing up, in this Letter we have proposed a method 
to create a generalization of long-range Cooper pairs us- 
ing bosonic atoms in an optical lattice. We suggest to 



prepare a Mott insulator of Bell pairs and dynamically 
melt this state into the superfluid regime. A spin inde- 
pendent interaction guarantees that pairing is preserved 
and that the system becomes a gas of long-range cor- 
related pairs, which spread over approximately half the 
lattice size. This mechanism works even if the process 
is not adiabatic or the initial state does not have trans- 
lational invariance. Most important, our proposal rep- 
resents a natural extension of current experiments with 
optical superlattices @, 0] • 
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